Mon. Not. R. Astron. Soc. 000, (2001) Printed 1 February 2008 (MN I^TeX style file vl.4) 



Separation of instrumental and astrophysical foregrounds 
for mapping CMB anisotropics 



J. Delabrouille^ G. Patanchon^ & E. Audit^ 

^PCC — College de France, 11, place Marcelin Berthelot, F-75231 Paris, France 

^ CEA/DSM/DAPNIA Service d'Astrophysique, CEA/Saclay 91191 Gif-sur-Yvette Cedex, France 



o 
o 

CO 
(N 



> 

00 

On 
O 

O 

o 



Accepted 2001 January 00. Received 2001 January 00; in original form 2001 January 00 



ABSTRACT 



For the most sensitive present and future experiments dedicated to Cosmic Mi- 
crowave Background (CMB) anisotropy observations, the identification and separation 
of signals coming from different sources is an important step in the data analysis. This 
problem of the restitution of signals from the observation of their mixture is classi- 
cally called "component separation" in CMB mapping. In this paper, we address the 
general problem of separating, for millimeter-wave sky mapping applications, compo- 
nents which include not only astrophysical emissions in two-dimensional maps, but 
also one-dimensional instrumental effects in the data streams. We show that compo- 
nent separation methods can be adapted to separate simultaneously both astrophysical 
emissions and components coming from time-dependent foreground signals originating 
from the instrument itself. Such general methods can be used for the optimal process- 
ing of low-redundancy observations where multi-channel observations are a precious 
tool to remove systematic effects, as may be the case for the Planck mission. 
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1 INTRODUCTION 

Mapping and interpreting sky emissions in the millimeter 
and sub- millimeter range, and in particular Cosmic Mi- 
crowave Background (CMB) temperature anisotropics, is 
one of the main objectives of present and upcoming obser- 
vational effort in millimeter- wave astronomy. This is, in par- 
ticular, th e objective of balloo n and space-b orne missions as 
Archeops ( Benoit et al. 2001), Boomer ang (|de Bernardis et 
al. 19991) , Ma xima (|Hanany et al. 200Q| ), TopHat (Bilverber^ 
200CI), MAP (IWright 1999|) and Plan"ck ([Lamarre et al. 2000[ 



these n 



required for precision tests of the cosmological models and 
measurements of the cosmological parameters (the Hubble 
constant ifo, the total matter density 11m, the cosmologi- 
cal constant A, . . . ), however, is such that it is necessary to 
disentangle in the data the contribution of several distinct 
astrophysical sources, all of which emit radiation in the fre- 
quency range used for CMB observations, i.e. the 10 GHz - 
800 GHz range. 

This problem, now traditionnally called "component 
separation" by the CMB commu nity, has been adressed 
already by a number of authors ( Tegmark Efstathiou 



Bersan^lli et al 20001). Among the scientific objectives of igggj. [Hobson et al. 1998|; iBouchet &: Gispert 1999|; iBouchet 



issions, tne precise measurement ot primordial tem- Pmnet Sethi 1999|; iBaccigalupi et al. 2QQQ| ; [snoussi et 



peratur e and/ or polarisation tluct nations ot tlie (JIVIB is one 2001 



of the priorities. 

The importance of measuring anisotropics of the Cos- 
mic Microwave Background (CMB) to constrain cosmolog- 
ical models is now well established. In the past ten years, 
tremendous theoretical activity, motivated by upcoming am- 
bitious balloon-borne and space missions, demonstrated that 
measuring the properties of these temperature anisotropics 
will constrain drastically the seeds for structure formation 
as well as the cosmological parameters describing the mat- 
t er content, the geometry, and the evolution of our Universe 
( Hu Sugiyama 1996| ; Jungman et al. 1996| ) . The accuracies 



Maino et al. 20Ql|) using different assumptions and 
analysis techniques. It is a particular application of the prob- 
lem of "source separation", also called in signal processing 
the "cocktail party" problem, in which a number of source 
signals (timestreams, images. . . ) has to be recovered from 
the observation of their mixture in a number of captors. 

So far however, only the separation of astrophysical 
components has been adressed using a component separa- 
tion formalism. In addition, the separation has been imple- 
mented for ideal white noise on the observed maps. Real 
measurements however often suffer from instrumental ef- 
fects, coloured noise, and involve the reprojection of one- 
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dimensional data streams onto somewhat irregularly sam- 
pled two-dimensional maps, where some pixels may be ob- 
served several times and/or others not covered at all by 
the observations, depending on the resolution of the recon- 
structed map. In this paper, we investigate the possibility 
to extend the W iener solution discussed by Te gmark & Ef- 
stathiou ( 1996| ) and Bouchet & Gispert (199£) to deal also 
with instrumental foregrounds. We test the method in more 
realistic situations which include non-white detector noise 
and uncertainties induced by the reprojection of timelines. 

The remaining of this paper is organised as follows: In 
section ^, we introduce the general formalism of component 
separation and discuss the possibilty to include instrumental 
components and their specificities. In section ^, we discuss 
the simulation of Planck observations containing the super- 



position of CMB and thermal dust emission, an instrumen- 
tal effect, low-trequency detector noise, and white noise. In 
section kl we discuss the Wiener method implemented to 
separate these components in the simulated Planck data, 
and show the results we obtain. In section ^ we discuss the 
results and the generality of the method. We conclude in the 
last section. 



2 COMPONENT SEPARATION 

One of the key aspects of the component separation is the 
exploitation of the presumably known (at least to a large 



extent) dependence of the component emission on the fre- 
quency ot observation v. Although it is not alwa ys necessary 
that the frequency scalings be known a priori (|3accigalupi 



irequ 

et al. 2QQq| ; pnoussi et al. 2QQl| ; [Maino et al. 2QQl| ), sepa- 



ration techniques rely on the existence of such a frequency 
scaling. In this framework, for multifrequency observations 
of the sky, the observed emission at frequency v is the linear 
superposition of the emission of several components: 



/(n, 



(1) 



where Ic{n,iy) stands for the emission of component c as a 
function of direction ft and frequency u. For each component, 
the emission can be well approximated as the product of a 
frequency spectrum Ac{iy) and a spatial template A(n). 



/c(n,zy) = Ac{iy)Ac{n) 



(2) 



This approximation is very good for small portions of the 
sky. For larger maps (especially covering both low and high 
galactic latitudes), there can be variations of the emission 
law of the galactic components. 

The decomposition in equation ^ is not unique, as a 
multiplicative constant factor can be taken from A into A. 
For convenience, we adopt the convention that Ac is the 200 
GHz template of component c in units of AT/Tcmb, that 



/c(n,200 GHz) 



^(TcMB,200 GHz) 



(3) 



TcMB 

with TcMB = 2.726 K, and Bi, is the Planck Blackbody 
spectrum. 

For the present application, we will assume that the 
frequency spectra of astrophysical component emissions are 
known, which is a reasonable assumption in the 100-800 GHz 



frequency range, where the sky emission at moderate galac- 
tic latitude is dominated by CMB anisotropics and thermal 
galactic dust emission (see figure The CMB anisotropy 
spectrum is the derivative of a blackbody with respect to 
temperature: 

AcMB (ly) (TcMB , ly) (4) 

Thermal dust emission is well fitted in this frequency range 
by a thermal spectrum with a emissivity and a tempera- 
ture of 17.5 K: 



(^) 



( ^ V 

V200 GHz; 



^.(17.5K,z/) 
B^(17.5K, 200 GHz) 



(5) 



This temperature of 17.5 K can cha nge slightly (at the level 
of a few Ke lvin) across the sky (Schlegel, Finkbeiner & 



Davies 1998), which is a motivation for a local treatment of 



component separation (as compared to full-sky methods). 

In addition to the known frequency dependence, one 
can introduce priors on the spatial statistical properties of 
the foregrounds. For Gaussian random fields, exploiting the 
prior knowledge of the autocorrelation leads to methods as 
the Wiener solution implemented and tested on simulations 
by Te gmark & Efstathiou (1996) and Bouchet & Gispert 
( 119991 ). 

In addition to diffuse galactic and extragalactic compo- 
nents, several populations of point-like sources are expected 
to contribute to the microwave sky emission. De dicated tech- 



nique s are being developed for their extraction ( Vielva et al 



2001). At this point however, techniques developed to sep- 



arate components in CMB experiments did not address the 
problem of non-astrophysical components, i.e. of signal fluc- 
tuations in detector readouts that are not due to astrophys- 
ical emission, fixed on the sky, but to the physical emission 
of other objects radiatively, or conduct ively, coupled to the 
detectors. 

One exemple of such a "foreground" is the thermal emis- 
sion of a slightly emissive telescope which temperature fluc- 
tuates in time. Detectors experience a time-dependent ad- 
ditive fluctuation of their signals, which, reprojected on the 
sky according to each detector's speciflc scanning strategy 
in the process of recombining timelines into maps, gener- 
ate fake position-dependent anisotropics on each detector's 
reconstructed sky map. 

Another example, speciflc to bolometric detectors, is the 
signal generated by temperature fluctuations of the thermal 
bath required to cool the detectors to sub-kelvin tempera- 
tures. Such fluctuations have been seen in detector timelines 



for the Archeops Balloon test flight (Benoit et al. 2001), and 
are expected to be present (at a much lower level) in Planck 
HFI data as well. 

Finally, one other interesting example is that of atmo- 
spheric emission for sub-orbital experiments, in which case 
time and space dependent fluctuations of the atmosphere 
temperature and composition in the line of sight generate 
fluctuations in the detectors timelines. 

Because of the process of their generation, such fore- 
grounds are seen by all detectors simultaneously. However, 
as flelds of views of all detectors of a CMB experiment as the 
current ones do not point towards the same point of the sky 
at the same time, they cannot be treated as another com- 
ponent on the sky. It is a component which "lives" in the 
time domain, and which is reprojected on the sky according 
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to the scanning strategy, in a detector- dependent way. As 
a consequence, a given pixel on the sky, being observed at 
different times by different detectors, is not pohuted by the 
same instrumental foreground for different detectors. 

Quite often, such instrumental foregrounds are essen- 
tially low frequency drifts, although one might imagine scan- 
synchronous temperature fluctuations depending on the ori- 
entation of the experiment with respect to external sources 
of heat as the Sun, or line-spectra fluctuations induced by 
cyclic power dissipation or by vibrations within the experi- 
mental set-up (e.g. engines). 

In the case of random slow drifts, a common way to 
handle this kind of systematics is simply to filter them out 
high paoD filter. Thio method hao the merit of oim — 



usmg a 

plicity pnd efficiency^ but reduces the sensitivity of the ex- 
periment to scales larger than ojscan x /filter, where ^scan IS 
the scanning speed on the sky and /filter the filter cutout fre- 
quency. For experiments aiming at measuring anisotropics at 
all available scales, including the largest ones, more sophis- 
ticated destriping methods are implemented. Such methods, 
whatever the exact algorithm used, essentially rely on redun- 
dancies (i.e. coming back to the same point of the sky on 
different timescales) and some prior assumptions about the 
slow drifts. They are not efficient in removing fluctuations of 
the systematic effect on timescales smaller than the typical 
"redundancy timescale" , unless they have built in some fil- 
ter which effectively does the equivalent of a high-pass filter 
on the timelines. 

For a scan strategy such as that of Planck, a 
timeline- proces sing destriping te chniqu e as that of Cot- 
tingham ( 1987 ) and Delabrouille ( 1998 ), which would non- 
destructively reproject "near-optimally" each single timeline 
onto a map of the sky while removing slow drifts of various 
origin, may be quite unefficient in removing stripes in re- 
gions near t he ecliptic plane, simply because of a lack of 
redundancy ( Delabrouille et al. 2000| ). The Planck scanning 
indeed is such that each detector's field of view scans the 
sky on large (~ 85 degree radius) circles centered near the 
ecliptic. The center of the circles is shifted by a few arcmin- 
utes every hour or so, the final field-of-view trajectory for 
each detector being a bundle of large circles intersecting each 
other near the ecliptic poles, with no intersection, or very 
few, in the region of the ecliptic (see figure 0). 

In addition, removing correlated low-frequency drifts in- 
duced by the detector pickup of slow thermal fiuctuations 
on the pay load would be non-optimal, because it does not 
take advantage of the physical origin of the systematic which 
gives rise to correlations. Taking care of this by usi ng map 
making techniques as adapted versions of MADCAP (Borrill 
1999) \^ hich would treat all detectors simultaneously, taking 
into account noise (i.e. systematic effect) correlations is a 
possibility, but is very intensive computationnaly. In this pa- 
per, we investigate the option which consists in treating the 
time-domain systematic effect as another component, "liv- 
ing" in the time domain and not on two-D maps, taking care 
of the specificities this implies as compared to sky- attached, 
more traditionnal, astrophysical foregrounds. 



3 SIMULATIONS 

We now investigate specifically the separation of instrumen- 
tal and astrophysical components on simulations tailored for 
the Planck mission. Much of the discussion can be straight- 
forwardly adaptated to other experiments as Archeops or 
Boomerang. The Planck mission, scheduled to be launched 
by ESA in 2007, will measure full-sky temperature and po- 
larisation anisotropics of the sky emission in the frequency 
range 30 to 850 GHz. The data will be taken using two in- 
struments, the Low Frequency Instrument (LFI), and the 
High Frequency Instrument (HFI). The LFI consists of 56 
radiometers cooled to 20K, split into four frequency channels 



in the 30-1 00 GHz range (Mandolesi et al. 2000; BersaneUi 



et al. 20QO|) , and the HFI of 48 b olometers split into 6 chan- 



nels between 100 and 850 GHz ( [Lamarre et al. 2000| ). This 
large frequency coverage, associated to unprecedented detec- 
tor sensitivity, will allow to measure microwave anisotropics 
with a final sensitivity of about 6 /xK per 8 arcminute pixel, 
provided all systematic effects remain under control. 

We will investigate a component separation method 
that could be used for this mission among others, restrict- 
ing ourselves to two sky components (CMB and thermal 
dust emission), one low frequency time-domain systematic 
effect due to the pickup of thermal emission from a part 
of the instrument, and realistic noise, including both white 
noise (in the time domain) and low-frequency detector noise 
(with a knee frequency of 0.01 Hz, as expected for Planck 
HFI detectors). We will process the observations of the 48 
HFI detectors only (6 channels) to avoid overconstraining 
excessively the problem. The method can be generalised to 
include LFI detectors as well, and more astrophysical (e.g. 
free- free emission, synchrotron, S unyaev-Zel'dovich (SZ) ef- 
fect dSunyaev Zel'dovich 197^ ), another dust component, 
zodiacal emission) or time-domain (e.g. temperature of ref- 
erence load, thermal bath fiuctuations) components. 

We implement the method on small square maps lo- 
cated near the ecliptic plane, convenient for a Fourier pro- 
cessing of the data. The motivation for the choice of location 
is that in regions near the ecliptic poles, redundancies in the 
measurements are such that standard "optimal" map mak- 
ing or destriping methods may be sufficient to remove the 
systematic effect, whereas we do not anticipate that this will 
be the case near the ecliptic. We pick a region located at co- 
ordinates a = 204 degrees and ^ = 11 degrees (A = 198 
degrees and p = 19.5 degrees), close to the ecliptic plane 
(for the scans to be almost parallel) and at a high galactic 
latitude (70 degrees). We choose, for the specific illustra- 
tion of our discussion here, a 272 x 272 pixels square map 
with two sides parallel to the equator, and a pixel size of 2.5 
arcminutes. 

The scanning strategy we adopt for generating our data 
is close to one of the possible Planck scanning strategies. It 
consists in spinning the spacecraft at 1 RPM so that the 
beam axis, for each detector, scans the sky along large cir- 
cles with 80 to 90 degrees angular radius (depending on the 
detector). The spin axis is held fixed for about an hour, so 
that each of the detector main beams scans the same cir- 
cle on the sky about 60 times. Every hour, the spin axis is 
displaced by about 2.5 arcminutes to follow the apparent 
yearly motion of the Sun, and get a complete sky coverage 
in about 6 months. Figure |l| illustrates the Planck scanning 
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North 
Ecliptic 
Pole 




South 
Ecliptic 
Pole 



Figure 1. Illustration of the Planck scanning strategy. On the left panel, the trajectory of the field of view of one single detector is 
shown. On the right panel, we show the trajectories of the centers of the fields of views of three distinct detectors for three spin axis 
position of Planck. For adjacent spin axis positions, the trajectories of all detectors are close to parallel for pointings near the ecliptic 
plane. 



strategy. Due to the detector setup in the focal plane, the 
scan- angle is slightly different for different detectors. 

We now turn to generating synthetic observations of 
this region for all Planck HFI detectors, featuring CMB, 
dust, an instrumental effect, and detector noise. 



3.1 Astrophysical components 

The CMB component is a COBE- normalised, randomly gen- 
erated realisation of CMB anisotropi cs obtained using 
predi cted by the CMBFAST software (Seljak & Zaldarriaga 
2000) for a Universe with Hq = 65 km/s/Mpc, Q.m = 0.3, 
= 0.045, A = 0.7. The CMB anisotropy emission is scaled 



as a function of frequency using the temperature derivative 
of the Planck law at T=2.726 K. 

The dust component is simply the extrapolation of the 
IRAS observed 100 micron map (bright stars removed) with 
an assumed emission proportional to x Bt^(T — 17. 5i^). 
The sum of the two for each frequency channel gives the sky 
emission observable in that band in our simulation. 

Input templates for CMB and dust emission can be seen 
in figure |5|. 

We observe the simulated maps assuming locally a par- 
allel scanning strategy such that the scans are distant, for a 
given detector, of 2.5 arcminutes from each other. We adop t 
the focal plane geometry described in ( Lamarre et al. 2000 ). 
For each detector, a purely astrophysical timeline is gen- 



erated for each Planck detector by observing the properly 
weighted superposition of CMB and dust maps in a set of 
samples obtained each by integrating the sky emission in a 
symmetric Gaussian beam, properly positioned with respect 
to the map according to our selected scanning strategy. The 
beam size ranges from 10 arcminutes (at 100 GHz) to 5 ar- 
cminutes (at 857 GHz). 

After this operation, we get a timeline in which only 
samples corresponding to observations in our generated map 
are non zero. At equal times, the 48 detectors do not point 
towards the same point on the sky because of the focal plane 
geometry. Shifts up to a few degrees, due to the angular 
extent of the image of the focal plane of the sky, separate 
the fields of views (FOVs) of the different detectors on the 
sky. 



3.2 Systematic effect 

A simulated instrumental foreground is generated in the 
form of a timeline corresponding to a plausible temperature 
fluctuation of an optical element: we generate a random, 
correlated gaussian stream, with a spectrum (as a function 
of frequency /) proportional to /~^. This is quite typical 
of temperature fluctuations observed in other experiments. 
We assume t hat a prior timeline treatment as, e.g., that of 
Delabrouille ( 1998 ) has permitted the removal of the effect 
of such fluctuations on timescales larger than the spinning 



© 2001 RAS, MNRAS 000, 



Separation of foregrounds for CMB mapping 5 



frequency of Planck, corresponding to one revolution of the 
satellite (one minute of time for a 1 RPM scanning). In the 
end, the simulated instrumental effect due to temperature 
fluctuations has a 1//^ spectrum between /spin 0.017 Hz 
and /sampling — 170 Hz, and falls off to negligible values 
outside this interval. 

The impact of the instrumental effect on each detector 
timeline depends on the value of a coupling coefficient for 
each detector. These coupling coefficients tells us how much 
the physical temperature fluctuations impact each detector 
reading. They depend on the emissivity properties of the 
object whose temperature fluctuates (as a function of fre- 
quency ly), and on the optical and geometrical properties of 
the setup. We assume here a z^^^^ dependence (close to the 
Ray leigh- Jeans limit of a greybody with emissivity propor- 
tionnal to as expected typically in our frequency range 
for aluminium) . The requirement on the temperature stabil- 
ity of the Planck pay load, including the mirrors, is 0.2 mK. 
Such a temperature stability can be achieved in the very 
quiet space environment of the L2 Lagrange-point orbit of 
Planck, but not very easily for balloon-borne experiments as 
Archeops or Boomerang. For our simulation, we assume a 
pessimistic temperature fluctuation of about 1 mK RMS in- 
tegrated on timescales between 0.017 Hz and 170 Hz (5 times 
higher than the expected upper limit for Planck), with the 
above 1//^ spectrum and with a 3 % total emissivity at 200 
GHz. 

This assumption of a strong instrumental effect, above 
expectations for Planck, make the dust contribution sub- 
dominant in all frequency channels. This allows us to test 
the separation of a "weak" astrophysical component. 

Note that no detector dependence of the coupling within 
a frequency band is assumed in our simulation. Such a de- 
pendence may exist in other cases as temperature fluctua- 
tions of a shield, for instance, where geometrical considera- 
tions have to be taken into account in the coupUngs. 

3.3 Detector noise 

For each detector, we generate a Gaussian coloured noise 
stream with a spectrum of the form Sn{f) = K x (1 + 
/knee//), i-G. a typical detector noise with a white part and 
a 1// part. The knee frequency (frequency at which the two 
parts contribute equally) is taken to be /knee = 0.01 Hz, 
which corresponds to laboratory measurements performed 
for the Planck HFI. The noise is assumed to be uncorre- 
lated between detectors. Noise flgures for the white noise 
level are representative of predictions for Planck HFI de- 
tectors. Again, we assume that very slow fluctuations have 
been removed, so that there remains no power outside the 
0.017-170 Hz frequency band. One noise timeline is gener- 
ated for each detector, and added to the appropriate data 
streams, which contain thus each the contributions from two 
astrophysical components, one instrumental component, and 
coloured noise. 



3.4 Reprojection on sky maps 

For convenience, we will perform our treatment on two- 
dimensional data rather than timelines. In the present im- 
plementation of the component separation, this allows to 



work easily in the Fourier space while taking into account 
the different detector pointings, in small regions of paral- 
lel scanning. We reconstruct a two-dimensional "observed" 
map for each detector by fltting to the observed timelines 
along with corresponding pointings a quintic surface defined 
in each pixel, matching best the observed data. A map of the 
systematic effect alone, obtained in the same way from the 
systematic effect timeline of a fictitious detector centered on 
the optical axis of Planck, is shown in figure ^ 

In figure |^, we display the reconstructed observed maps 
for one single detector of each channel. These maps include 
the reprojected noise and systematic effect for each detector, 
as well as the observed astrophysical components. In figure 
^ we display the observations for 3 detectors of the 217 
GHz channel, which show the detector dependence of the 
observed systematic effect which, being synchronous on the 
timelines, is not coincident on the maps because of the shifts 
between detectors. 

In table 0, we give for all channels the ratio between the 
RMS of the signal (CMB or dust) and the error performed 
on its measurement with our simulation setup. The error 
contains contributions from detector noise, from the other 
astrophysical component, from the instrumental effect, and 
from observational errors (beam smoothing and reprojection 
effects) . 



4 IMPLEMENTATION OF A COMPONENT 
SEPARATION 

The separation of astrophysical components we implement 
is possible essentially because the components have distinct 
emission spectra as a function of radiation wavelength. This 
is illustrated in figure ^. If the number of components and 
their respective spectra were perfectly well known, and if 
measurements were noiseless, multifrequency observations 
(with at least as many observations as there are compo- 
nents) would allow in principle the perfect recovery of the 
emission of each of the component at each point. In practice, 
imperfections of the measurements and the model make the 
problem slightly more complicated. 

4.1 First principles 

In the case where astrophysical components only are present, 
we are given a number of observations of the same sky 
(or fraction of sky) at different frequencies z/, each of which 
can be written as a linear superposition of emissions from 
the various astrophysical sources of radiation (CMB, dust, 
synchrotron, SZ. . . ) with some noise added: 

(x) = A^cTc (x) + (x) (6) 

This equation can equivalent ly be written in Fourier space, 
as 

(^) = A^cTc (k)+n^{k) (7) 

Auc is the "mixing matrix" of this component separation 
problem, which is determined by the spectrum of the emis- 
sion of the astrophysical components (up to a global calibra- 
tion factor for each component) and by the optical response 
of the instrument. 
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353 GHz 545 GHz 857 GHz 



Figure 2. Observations for one detector at each wavelength. For illustrative purposes, we have chosen here a strong systematic effect, 
which is comparable to CMB emission at low frequencies, and dominates over dust at high frequencies. 




Figure 3. Observations by three different detectors of the 217 GHz channel. In our example here, the observations at this frequency 
are dominated by CMB anisotropics and stripes induced by the instrumental effect . The detector corresponding to the middle panel is 
polarised and thus slightly noisier after recalibration of the observations, because half of the input power (here assumed to be unpolarised) 
is rejected. The increased granularity due to high frequency noise can (barely) be seen on the observation. The detector-dependent shift 
of the systematic effect with respect to sky emission is clearly visible on these observed maps (look at the distance of the darkest stripe 
to the upper left corner of the maps). 
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Table 1. Signal to Noise ratio for CMB and dust for observations in all channels. Here we consider as "noise" the 
error from all sources : detector noise, smoothing with the beam, other components, and reprojection effects. 



channel 100 143 143p 217 217p 353 545p 857 



number of det. 


4 


3 


9 


4 


8 


6 8 6 


S/N CMB 


1.4 


1.2 


0.95 


0.77 


0.65 


0.16 0.01 10-^ 



S/N dust 8.9 X 10-3 2.0 x 10"^ 1.7 x 10"^ 5.9 x 10"^ 5.2 x 10"^ 0.14 0.22 0.26 



.-3 



.-5 



Synch rotron 



CMB 



Thermal SZ 



Systematic 

Spinning dust (?)\ ' 



10 100 

Frequency (GHz) 



1000 



Figure 4. Foreground emission in the frequency range of inter- 
est, as compared to CMB anisotropy spectrum. The existence of 
three major galactic foregrounds, which are thermal dust emis- 
sion, synchrotron, and bremstrahlung (also called free-free, not 
represented here), is well established. There is some evidence to- 
wards the existence of microwave emission due to rotating dust. 
SZ effect generates foreground emission towards clusters of galax- 
ies, and possibly filamentary structure between them. We have 
added a plausible spectrum for a systematic effect due to tem- 
perature fluctuations of an aluminium telescope mirror. Relative 
levels are arbitrary, so that the curve for each foreground emission 
can move up and down depending on the level of contamination. 
The assumed shape of a p ossible spinning dust component is from 
Draine & Lazarian dl998t) 



If there are more observing frequencies than astrophys- 
ical components, and if the distributions for n and T are 
Gaussian and their autocorrelations N and S are known, 
the Wiener solution can be written as: 

f = [S~^ + A'N-^A]-^A'N-^ m = Wm (8) 

Equation ^ defines the Wiener matrix W — [S~^ + 
A^N~^A]~^A^N~^, which gives the best estimate f of the 
sky temperature T. 

4.2 Including instrumental effects 

Let us now assume that a number of instrumental effects 
are to be added to our model. The measurement for each 
detector can be written as: 



TTid {x) = AdcTc {x) + A'dzSz (x, d) + Ud {x) 



(9) 



where now d indexes detectors. Sz{x,d) is the detector- 
dependent projection of the instrumental component. This 
"reprojection" of the time- dependent instrumental compo- 



nent is not the same for all detectors because although all 
detectors experience the same instrumental effect (up to 
detector-dependent coupling factors) at the same time, they 
do not see the same pixel at the same time. Therefore, the 
systematic effect which pollutes the measured map in x is 
not the same for two different detectors d and d' . Because of 
this, it is not possible to recombine simply the measurements 
of all the detectors at a given frequency into one single map. 

In all generality, Sz (x, d) depends on the scanning strat- 
egy and on the relative orientation of the focal plane and 
the scanning direction. Therefore the relationship between 
reprojected effects Sz{x, d) and Sz{x^ d') between two differ- 
ent detectors d and d' may be far from trivial (see appendix) . 
However, if we restrict ourself to region of the sky where all 
the scans are parallel, there is a constant shift Ax^ between 
the direction pointed by detector d and that of the optical 
axis, which is set entirely by the geometry of the image of 
the focal plane on the sky. In that case equation (^ can be 
rewriten: 



md{x) - AdcTc{x) + A'dzSz{x - Axd) + nd{x) 



(10) 



where Sz{x) is now the instrumental noise projected on the 
sky in the direction of the telescope optical axis. Although 
it is not possible to average in bands without loosing infor- 
mation on the systematic, the system has a linear structure 
close to that of equation ^, except for the additionnal sys- 
tematic effect. 

Rewriting equation ^ in Fourier space, we get 

rridik) = AdcTcik) + A'dz e^p{ik.Axd)Sz{k) + nd{k) (11) 

This has exactly the same structure as that of equa- 
tion ^ if we include the spatial frequency dependent term 
exp(iA;.Axcz) into matrix A\ and regroup matrices A and A' 
into one single matrix (and {Tc,Sz) into one single vector). 
Therefore, we can apply a Wiener method to our extended 
problem. There are a few technical dif ferenc es with the im- 
plementation of Bouchet and Gispert (1999): 



• in the present case, because of the phase shift for the 
systematic component, the mixing matrix A is complex; 

• because of the problem asymmetry induced by the priv- 
ileged scanning direction, the Wiener matrix W depends on 
k not simply on |^|; 

• the noise we assume is not white. 

• we have one observed map for each detector, not one 
single map per frequency band. 

• the observed maps are not simply the superposition of 
the true pixellised maps with white noise added, but are 
reconstructed from timelines and pointings by a global fit to 
the data, which takes into account some of the inaccuracies 
of the map-making. 
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4.3 Finite beam effect 

The effect of symmetric finite beams is to convolve the maps 
with a symmetric kernel B(x). For asymmetric beams, the 
integration of the sky in the beam is scan-direction depen- 
dent, so that the operation is not strictly a convolution, ex- 
cept in the case of local parallel scanning. In Fourier space, 
the sky temperature Tc is multiplied by the Fourier trans- 
form of the asymmetric beam pattern (properly oriented). 
The effect of the beam can therefore be taken straightfor- 
wardly into account by putting its ^-dependent effect into 
matrix A (but not matrix A\ as the systematic effect is not 
affected by the shape of the beam) . Each component Adc of 
matrix A needs to be replaced in equationjjj by the prod- 
uct AdcBd{k) (no summation), where Bd{k) is the Fourier 
transform of the beam map for detector d. We include the 
effect of the symmetric Gaussian beams in the evaluation of 
the mixing matrix A for the inversion in Fourier space. 



4.4 Implementation 

The implementation of the Wiener method in Fourier space 
requires the apodisation of the maps to avoid border effects. 
We take an apodisation function which goes to zero over 
15 pixels (out of 272) on each side of the map, following a 
sine function from 7r/2 to 0. This apodisation induces little 
smoothing on the spectrum in /c- space, and leaves the inner 
region of the map untouched. 

The implementation of the Wiener solution also requires 
prior information of the spectrum in /c-space, S{k)^ for all 
components. For the CMB, we assume the true underlying 
spectrum (which differs slightly from the actual spectrum 
of our realisation because of the sample variance). For dust, 
we assume a 1/(|£| + C)^ law, where C is a constant chosen 
to fit best the spectrum on the IRAS map. For the system- 
atic, we take the average of 10 realisations of a systematic 
according to the same generation law that we use for gen- 
erating the actual observed systematic. Therefore, for all 
components, we assume a spectrum which is close to the re- 
ality, but not exact (either because of sample variance, or 
because of different realisation, or a prior empirical fit). We 
consider that this is representative of the uncertainty we will 
have when dealing with actual measurements. We discuss in 
more details the effect of uncertainties in the model of the 
measurement (and in particular uncertainties in pointing, 
frequency scalings) in the next section. 

The simulation and inversion we implemented takes a 
few minutes of running time for 300 x 300 2.5 arcminute 
pixels maps, on a very modest workstation. 



4.5 Results 

In figure ^ we display the result of our inversion. Apart for 
border effects which here cannot be avoided, the original 
maps are well recovered. Border effects are due both to the 
apodisation and to discontinuity effects of the instrumen- 
tal effect map. Apodisation effects appear only in a Fourier 
implementation, and could be avoided in real space imple- 
mentations, at the cost of enormous required computation 
time. The discontinuity effect on the other hand is not a 
problem of the implementation. It appears because the map 



Table 2. Final Signal to Noise ratio for CMB and dust, for both 
the 48-detector and the 6-detector inversions. Figures in paren- 
theses are obtained when the output is compared with the original 
map smoothed with the effective beam. 

component CMB dust 

48-detector separation 7.0 (10.1) 3.5 (5.1) 

6-detector separation 4.1 (5.5) 2.5 (3.7) 



of the instrumental effect is not periodic, and the actual in- 
strumental effect limited to the region of observation is thus 
not exactly the same for all detectors, because of detector 
shifts. The effect remains small as long as detector shifts are 
small compared to the size of the map. 

The RMS error on CMB anisotropics measurement, 
computed as the difference between the true input CMB 
maps and the output map (obtained from the reprojection 
of timelines and suffering from convolution with the beam, 
additive noise, and possibly residuals of other components 
after the separation process) is AT/T 5.3 x 10~^ in the 
central region of the map, corresponding to a signal to noise 
ratio of about 7. 

Dust emission is very subdominant in all channels of our 
simulation, because we have chosen quite a strong instru- 
mental effect. This is unlikely to be the case for Planck, but 
may be typical of some balloon-borne experiments. Despite 
of the fact dust can barely be seen in individual detector 
observations, its emission, in the central region, is recovered 
with a signal to noise ratio of about 3.5, and a AT/T of 
5.3 X 10"^ (at 200 GHz). 

For both components, final results are much better than 
what could be expected from simple weighted averaging of 
all measured maps, as is clear when comparing the final S/N 
rates with figures from table |l|, and output maps (fig. ^ with 
maps of the observations |2| 

The stripes due to the instrumental effect are well re- 
covered, with a signal to noise of 7.2 in the central region of 
the map. Of course, the sides of the map are not recovered 
as well, because of both the apodisation and the non-toric 
geometry of our maps, which invalidates the Fourier phase 
shift approximation of the actual dispacement of each de- 
tector at the borders. 

In table |2| we show the signal to noise ratio obtained at 
the end of the inversion, both for the 48-detector inversion 
and for a 6-detector inversion with one detector per channel 
only. For each case, we give S/N computed in two ways. 
The first way is to compute the error from the difference 
between the "true" input map and the recovered map. The 
second is to compute the error from the difference between 
the input map smoothed with the effective beam (see |4.6[ ), 
and the recovered map. Values obtained with the second 
method are given in parentheses. The difference between 
the two is due to the part of the true astrophysical signal 
which is not recovered because of the imperfect resolution. 

The difference between the signal to noise ratio for the 
6-detector inversion and the 48-detector inversion is com- 
patible with the loss of sensitivity due to the reduction of 
the number of detectors by a factor of 8. 

In summary, the separation algorithm we implemented 
performs very satisfactorily, and can be a useful tool to con- 
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Figure 5. Results of the inversion. The input maps are displayed in the first column. Maps for CMB and dust are original maps, whereas 
the map for the instrumental effect is a reprojection of the effect's timeline. The output of the inversion, to be compared with the original 
maps, are displayed in the second column. The difference is shown in column three. For each component (on each line), the three maps 
are displayed with the same greyscale table (shown on the right of the figure) where figures give values for AT/Tcmb at 200 GHz. 
Signal to noise ratios of the final measurement in the center region (delimited by black lines on the difference maps are given below the 
difference maps.) 



strain instrumental effects and separate various components 
in CMB data for a scanning experiment such as Planck. 



4.6 Effective beams 



As T = Wm, we have T WAT + (noise term), W being 
the Wiener matrix and A the mixing matrix. This permits to 
compute an "effective beam" for each component, which is 
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Table 3. Effective co-scan and cross-scan beam sizes for CMB 
and dust, for both the 48-detector and the 6-detector inversions. 



Direction 



cross-scan 



48-det. CMB effective beam 7.99 arcmin 7.99 arcmin 
48-det. dust effective beam 6.83 arcmin 7.07 arcmin 



6-det. CMB effective beam 
6-det. dust effective beam 



10.0 arcmin 10.0 arcmin 
9.82 arcmin 12.9 arcmin 



the backward Fourier transform of the appropriate diagonal 
term oiW{k)A{k). 

The final effective beam is due to a combination of the 
instrumental beam (which gets into matrix A) and of the 
attenuation of some spatial frequencies, where noise domi- 
nates, by the Wiener filter. 

We compute effective beams for CMB and dust for both 
the 48-detector and the 6-detector inversions. These effective 
beams are well fitted by two-dimensional elliptical Gaus- 
sians. Because of the stripes, which increase the effective 
noise term on high frequencies perpendicular to the scan- 
ning, but not parallel to the scanning, beams are slightly 
elongated perpendicularly to the scanning, but this is not 
significant for the CMB. It is quite significant for the dust 
(especially in the case of the six-detectors inversion). This 
was to be expected, since the dust emission has a frequency 
dependence close to that of the systematic effect and a much 
lower amplitude. Effective beam width for both inversions 
in co-scan and cross-scan directions are shown in table |^. 

The effective resolution depends on the number of de- 
tectors because in the case of 6 detectors, high frequency 
noise starts to dominate the signal at lower k values than 
when 48 detectors are present, which leads to effective beam 
degradation because of the Wiener filtering, which weights 
modes with a S/{S + N) factor. Results shown in this table 
are pessimistic for Planck, and especially so for dust because 
the instrumental effect, which dominates at high frequencies 
where dust is the highest astrophysical contribution, is quite 
high as compared to what is expected. 



5 DISCUSSION 

Our simulations have shown that parallel stripes due to the 
simulated systematic effect can be separated very well from 
astrophysical emission in the particular case we studied here. 
We now discuss the reasons why it works and the dependence 
of the result on specific assumptions. 



5.1 Impact of the number of components to be 
separated 

The test of our method as described in this paper has been 
done in a particularly favorable situation as far as the num- 
ber of components is concerned. We have 6 channels and 48 
detectors for 3 components only. The system is thus widely 
overconstrained. There is, however, no fundamental prob- 
lem to adding more astrophysical com ponents, as already 
demo nstrated by Gispert and Bouchet ( Bouchet Gispert 
199£). Here, we simply demonstrate the possibility to treat 



instrumental effects as additionnal components, provided 



the inversion is made on a system encompassing all detec- 
tors, and not co-averaged frequency channels. 

We have performed tests of the method with the same 
number of components (two astrophysical and one system- 
atic) but with a very limited number of detectors and chan- 
nels (as little as three frequency channels and six detectors), 
with satisfactory results. 



5.2 Impact of the spectral properties of the 
systematic effect 

The spectral properties of the astrophysical components and 
of the noise is an important issue regarding the separation 
procedure. 

Let us imagine, for instance, that the systematic effect 
has almost the same frequency dependence as dust emission, 
e.g., an emission proportionnal to . This may happen for 
atmospheric emission in some cases of ballon-borne experi- 
ments, or in the case of Planck if some element of the space- 
craft, which temperature fluctuates and which is coupled 
optically to the detectors, has a greybody emission with 
emissivity. In this case, the frequency dependence does not 
help separating the two, and if all detectors were looking at 
the same point of the sky at the same time, there would be 
a very strong degeneracy between the dust and the thermal 
emission. The best we could do in that case would be to 
filter stripes out of dust maps only based on spatial prop- 
erties or on statistics (loosing at the same time some high 
frequency in the information on dust emission, in direction 
perpendicular to the stripes). However, if as in the method 
we propose, the information concerning the different point- 
ing direction of each detector is used, this degeneracy can 
be lifted to some extent. 



5.3 Further developpements 

The results we have presented here used a symmetric beam. 
However, it is straightforward to extend our procedure to 
an asymmetric beam. As long as scans are parallel, the ef- 
fect on the maps of an asymmetric beam can be written, 
mathematically, as a convolution. The Fourier components 
of the sky emissions are multiplied by the Fourier transform 
of the beam (properly oriented with respect to the scanning 
direction), which modifies their spectra. This can be taken 
into account easily in the formalism of the separation, where 
the priors on components depend on k and not simply on 

lfc|. 

The generalisation of our method to non-parallel scans 
is more difficult. Our component separation method as im- 
plemented here works only for parallel equidistant scans. 
In principle, small deviations from this assumption are not 
critical, as long as the deviations remain small as compared 
to the beam sizes. The potential problem when the scan- 
ning departs from this ideal case is not the distortion of 
the assumed spectra for the Wiener solution, but the ac- 
curacy of the model of a detector-dependent constant dis- 
placement which permits to write the system in the simple 
form of equation |ll|. We are currently working on relaxing 
the parallel-scanning assumption. 

A first possible generalisation concerns a portion of the 
sky where two sets of parallel scans cross each other (which 
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is often the case in actual CMB experiment). In this case, it 
is possible to adapt our formalism by adding two systematic 
effect maps (one for each scanning direction) instead of one. 



6 CONCLUSIONS 

In this paper, we have investigated the possibility to use 
Wiener filtering to separate astrophysical components and 
systematic effects producing slow drifts in scanning CMB 
anisotropy experiments as the Planck mission. We have 
shown that the "component separation" formalism applied 
so far only to astrophysical emissions can be generalised to 
take into account and separate instrumental effects as well. 

The algorithm we implemented is especially useful to 
remove systematic effects of known spectrum when there 
is little redundancy in the measurement for each detector 
channel, but high correlation between the systematic effect 
impact in the different detector channels. Whereas we have 
studied here the impact of one single systematic effect, that 
of slow temperature drifts of an optical element of the Planck 
telescope, the approach can be generalised and adapted to 
treat more effects, as for example one or more components 
due to atmospheric emission. 

One of the difficulties for the application of our method 
is to get a good evaluation of the elements of the "mix- 
ing matrix", i.e. matrix A, for the systematic effects. For 
Planck, an accurate modeling associated to numerical sim- 
ulations may permit to constrain these elements quite well, 
which may not be the case for other experiments/missions. 
Methods permitting to estimate simultaneously mixing ma- 
trix elements and the compon e nts themselves (Baccigalupi 
et al. 2QQQ| [Snoussi et al. 2QQl| |Maino et al. 2Q0l[ ) can have 
a natural application in the context we discussed here. 
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APPENDIX A: GENERAL FORMULATION OF 
THE FULL SEPARATION PROBLEM 

In this appendix, we write a general formalism for the sepa- 
ration of both sky-domain (astrophysical) and time-domain 
(instrumental) components for a scanning experiment. 

In a scanning experiment such as those currently used 
for CMB mapping, it is usual to write the data stream of a 
given detector d as a function of time as: 



St 



Mtpip + Tit 



(Al) 



Where st is the data stream for that detector (and t 
indexes time), Ip the brightness of a pixellised sky map in 
pixel p, and Mtp the pointing matrix, which tells us how 
much pixel p contributes to the data stream at time t. A 
summation convention is assumed over the repeated index 
p. nt is a noise stream, representing random noise, which 
may in some cases be coloured. 

From a component separation point of view, the sky 
brightness is the sum of several astrophysical components, 
Ip = AcTcp. Here, Ac is a single column of the "mixing ma- 
trix" of the component separation problem, corresponding 
to one single detector d. The full mixing matrix for all de- 
tectors is noted Adc- 

Putting now explicitely in the equation the detector de- 
pendence as an additionnal index and writing explicitely 
the summations to avoid confusions, we get: 



Sdt 



MdtpAdcTcp + Udt 



(A2) 



For a satellite experiment as Planck, the size of matrix 
Adc is about 100 x 6, and the size of matrix Mdtp about 
100 X 10^ X 10^, making it impractical to handle all the 
problems at the same time. As long as the noise of all detec- 
tors are independent, it is convenient to reproject timelines 
onto maps (one per detector, or possibly one per frequency 
channel) and then separate astrophysical components. 

Instrumental effects as the one we discuss in this paper 
can be added to the formalism, in which case formula |A2| 
becomes: 



Sdt = ^ ^ MdtpAdcTcp + ^ KdiXit + Udt 



(A3) 



In this formula, Mdtp is the pointing matrix for detector 
Adc the astrophysical component mixing matrix (which de- 
pends on the detector d essentially through the correspond- 
ing frequency channel), Tcp the sky emission of component 
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c in pixel p. xu are the instrumental components (tempera- 
ture fluctuations of spacecraft elements, atmospheric emis- 
sion. . . ) where i indexes components and t indexes time, 
and Kdi the mixing matrix for these instrumental compo- 
nents, which depends on the physical couplings between the 

sources and the detectors. Udt is the noise term. 

In principle, equation A3 (just as well as equation A2) 
can be rewritten in the general form s = AT + n. The form 
of equation A3, however, shows more clearly the physics 
underlying the model. 
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